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Abstract — A two-dimensional signal space for intensity- 
modulated channels is presented. Modulation formats using this 
signal space are designed to maximize the minimum distance 
between signal points while satisfying average and peak power 
constraints. The uncoded, high-signal-to-noise ratio, power and 
spectral efficiencies are compared to those of the best known 
formats. The new formats are simpler than existing subcarrier 
formats, and are superior if the bandwidth is measured as 90% 
in-band power. Existing subcarrier formats are better if the 
bandwidth is measured as 99% in-band power. 

Index Terms — Direct detection, intensity modulation, nonco- 
herent communications, power efficiency, spectral efficiency. 

I. Introduction 

INTENSITY modulation with direct detection (IM/DD) 
is widespread for low-cost optical communication sys- 
tems, e.g., wireless optical links 1T|-|I3| and short-haul fiber 
Unks im. IM/DD permits only the intensity of light to carry 
information. In contrast, coherent optical systems such as long- 
haul fiber links let data modulate the optical carrier's amplitude 
and phase via, e.g., Af-ary quadrature amplitude modulation 
(Af-QAM). Designing IM/DD formats with good power and 
spectral characteristics is challenging lH], |]2], 0-117]. 

In the absence of optical amplification, IM/DD systems 
can be modeled as additive white Gaussian noise (AWGN) 
channels with nonnegative inputs |1 Ch. 5], 121, Q, IS] 
Sec. 11.2.3]. Nonnegative Af-ary pulse amplitude modulation 
(Af-PAM) such as on-off keying (OOK) |1, Eq. (5.8)] is 
a natural modulation format but it is power inefficient for 
M > 2 |9|. Subcarrier modulation (SCM) allows using M- 
QAM by adding a direct current (DC) bias to the electrical 
signal to make it nonnegative [1 Ch. 5]. The DC bias does 
not carry information if it is independent of the transmitted 
information. A signal space for IM/DD channels was presented 
in Q and power-efficient subcarrier modulation formats were 
designed. In our prior work, a three-dimensional signal space 
for IM/DD, whose signal sets are denoted as raised-QAM \5i, 
was used to numerically optimize modulation formats for 
different power constraints ||6l, ||7|. 
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Fig. 1. Passband transceiver of IM/DD systems. 



In this work, we present a two-dimensional signal space 
for optical IM/DD systems. The resulting modulation formats 
have simpler modulator and demodulator structures than the 
three-dimensional formats studied before. Their power and 
spectral efficiencies are evaluated and compared to the pre- 
viously best known formats. 

II. System Model 

Consider an AWGN channel whose input x{t) is nonneg- 
ative. The symbols u{k), for k = . . . , — 1, 0, 1, . . ., are inde- 
pendent and uniformly distributed over {0, 1, . . . , M— 1}. The 
modulator maps each u{k) to a real and nonnegative waveform 
belonging to the signal set S — {so(i), si(t), . . . , SM-i{t)}, 
where s,{t) = for i ^ [0, T), i = 0, 1, . . . , Af - 1, and T is 
the symbol period. The transmitted waveform is 



x{t) = ^ Su{k)(t-kT). 

k— — oo 

The received signal is modeled as 

y{t)^x{t)+n{t), 



(1) 



(2) 



where n{t) is a zero-mean white Gaussian process with 
double-sided power spectral density Nq/2. The demodulator 
is a correlator or matched filter receiver with a minimum- 
distance detector, i.e., it minimizes the symbol error rate at a 
given signal-to-noise (SNR) ratio 1,10. Sec. 4.1] and puts out 
u{k) as the estimate of u{k). 

The AWGN model is reasonable for IM/DD systems when 
the dominating noise is from the receiver itself, and not from 
optical amplifiers [1 Ch. 5], H, Q, iS Sec. 11.2.3]. A 
passband model for IM/DD systems is depicted in Fig. [T] 
The electrical nonnegative waveform x{t) modulates a light 
source such as a laser diode. The information is carried by 
the intensity of the passband signal 



j(i) == ^/2cx{t) cos{2tt fot + 0{t)), 



(3) 



where c represents the electro-optical conversion factor in 
watts per ampere (W/A) [11], /o is the optical carrier fre- 
quency, and 9{t) is a random phase, uniformly distributed 
in [0, 2tt) and varying slowly with t relative to the symbol 
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Fig. 2. Two-dimensional constellations optimized for average and peak optical power Circles which are different between constellations of the same size 
are shaded differently. 



rate. The optical signal propagates through the optical medium 
depicted as an optical fiber in Fig.[Tl but which could also be a 
free-space optical link. The photodetector at the receiver, with 
responsivity r in AAV, detects the intensity of z{t). Without 
loss of generality, c and r are normalized so that the received 
electrical signal y{t) can be written as ^ B121 p. 155]. 

III. Signal Space Analysis 



The signals in S can be represented as 



N 



Si{t) = y^^ Sj^k (t>k{t) 



(4) 



fc=i 



for i = 0, ...,A/ - 1, where {(/-feWlf^i, N < M, is 
a set of orthonormal basis functions |5|. The vector repre- 
sentation of Si{t) with respect to these basis functions is 
Si = (si,i, Si, 2, ■ • ■ , Si,Af)- We may thus alternatively represent 
the signal set as S* = {sq, Si, . . . , Sji/_i}. 

Consider a two-dimensional signal space for IM/DD 
spanned by the basis functions 



where 



Mt) = 
Mt) = 

rect(i) 



1 ft 

— rcct — 
T \T 



2 ft 

- cos (27r/i) rect ( - 



1, ifO<i<l, 
0, otherwise. 



(5) 
(6) 



and the electrical subcarrier frequency / = 1/2T. In |i51 and 
our prior work ||6], ||7], a three-dimensional signal space with 
signal sets called raised-QAM was used to design modulation 
formats. This three-dimensional signal space is spanned by 
4>i (t) in (|5]l, and the in-phase and quadrature phase modulation 
formats' basis functions with / = l/T 17, Eqs. (13)-(14)]. 
The basis function 0i(t) represents the DC bias, and is used 
as in [51 to guarantee signal nonnegativity. 

We follow the same steps as in Theorem 1 in |7| for 
raised-QAM. The admissible region T, defined as the set 
of two-dimensional signal vectors satisfying a nonnegativity 
constraint, is a two-dimensional cone with vertex at the origin, 
an apex angle of cos~^(l/3) = 70.528°, and an opening in 
the dimension spanned by (/)i{t). 



A. Example Modulation Fonnats 

Fig- El presents several modulation formats designed for 
the admissible region T. For a certain power constraint, the 
formats are optimized to maximize the minimum distance 
between points. As in |5J, |7J, the average optical power 



Pr, 



lini — 
T^oo 2T 



z2(i)di 



-T 



lini — 

T^oo 2T 



x{t) di, (7) 



-T 



and peak optical power 



Po = max ■ 



2 t ^ ' 



(8) 



are used as design criteria. The modulation formats numer- 
ically optimized for average optical power are denoted as 
,'7p jv/, and for peak optical power as 3^ j^j. They are 
listed in Appendix |A] For illustration. Fig. |2] depicts the 
constellation points, regarded as the centers of unit-diameter 
circles, together with the admissible region T. 3^ is a 4- 
ary constellation optimized for both power measures. Together 
with ^/'p 3 and 3'p g, they are subsets of a lattice where the 
angle between its two basis vectors is cos^^(l/3), which is 
the apex angle of the cone. 

B. Performance Measures 

To evaluate performance, the uncoded and asymptotic (high- 
SNR) power gains with respect to OOK are considered [IJ, 
10, Q- The average optical power gain Po,gain over OOK 
for the same error rate performance is defined in [7. Eq. (32)] , 
and the peak optical power gain Po,gain with respect to OOK 
is defined in [7, Eq. (33)]. The average optical power is an 
important figure of merit for skin- and eye-safety in wireless 
optical links [1, Ch. 5], IJS], ||5], and the peak power measures 
tolerance against nonlinearities llT3l . 

The spectral efficiency measures the bit rate achieved in a 
given bandwidth. It is defined as 



77=^ [bit/s/Hz], 



(9) 



where Ri, = Rg logj M, Rg — l/T is the symbol rate in 
symbols per second, and W is the baseband bandwidth of 
x{t). In m, 161, Q, W was defined as the first null in the 
spectrum, i.e., the width of the main lobe, since most of the 
energy of a signal is contained in this main lobe. However, 
some modulation formats designed using the signal space in 
Sec. Hn] have no spectral nulls. Fig. |2] shows five formats that 



%,3 ^3 «i 



OOOK 



^P„.8 * 



I .5-1 



>QPSK 



"^PoAe * 



0^A,8 



*-P„,l6 #'^„,8. a^_ 



OOOK 



-1 - ♦ %,3 

-1.5- [>QPSK 



1^4 



(/^K 



* ^Po:8 



'«^A.,ifi 



.%.«* ♦■^^". 



* ®P„,16 



1.5 2 2.5 3 3.5 4 4.5 

7) [bit/s/Hz] 

Fig. 3. Asymptotic average optical power gain of modulation formats over 
OOK (K = 0.9). 
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Fig. 4. Asymptotic peak optical power gain of modulation formats over OOK 
(K = 0.9). 
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Fig. 5. Asymptotic average optical power gain of modulation formats over 
OOK (K = 0.99). 

lack spectral nulls. This makes this definition of bandwidth 
misleading. Instead, as in [3, p. 49], we will use the fractional 
power bandwidth W defined as the length of the smallest 
frequency interval carrying a certain fraction of the total power. 
Note, however, that [3, Eq. (3.15)] accounts only for the 
continuous spectrum of x{t). We will include both the discrete 
spectrum and the continuous spectrum. For IM/DD channels, 
the discrete spectral component at / = (at DC) represents 
the average optical power of a constellation 13, p. 47]. The 
fractional power bandwidth W is the solution to 
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Fig. 6. Asymptotic peak optical power gain of modulation formats over OOK 
(K = 0.99). 
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(10) 



where S^if) is the power spectral density of x{t), and K E 
(0,1). For constellations with uniform probability distribution, 
Sxif) can be obtained using |14, Eq. (3.7.6)], which depends 
only on the Fourier transform of the signals in S. 



IV. Performance Analysis 

Figs.[3]-(6]depict the (uncoded, high-SNR) average and peak 
optical power gains of our modulation formats with respect to 
OOK, and as a function of spectral efficiency 77. The fractional 
power bandwidth W was computed using K = 0.9 and 
K = 0.99. These choices of K are somewhat arbitrary but 
are commonly used. 

In addition to the modulation formats introduced in this 
paper, we consider subcarrier QPSK, nonnegative 4-PAM, 
and three-dimensional modulation formats from previous work 
optimized for average optical power C^/fp f^j) and peak optical 
power C^p j^) [7J. We next discuss the plots in Figs. [3]-(6] 

1) Po,gain VS. Tj with K = 0.9 (Fig. [i]).- For a fixed 
AI, modulation formats optimized for average power have a 
larger Po,gain than those optimized for peak power The three- 
dimensional constellation optimized for both power measures, 
■^4, has a similar r] and Po.gain as ."^p 3, whereas ,%^ has 



better power and spectral efficiency than '^p g. The three- 
dimensional constellations optimized for peak power, ^Sp g 
and "^p j^g, have higher r\ than ^€p^ g and '€p^ j^g. This comes 
at the price of a lower power efficiency. 3'p g has a higher r\ 
than S'p g, and 4-PAM has the highest r\ and lowest Po.gain- 
For K = 0.9, the modulation formats designed using the 
two-dimensional signal space have better spectral-efficiency 
characteristics than formats designed using raised-QAM. For 
example, S'p g has a better r\ than '£p -^g. 

2j Po.gatn'vs. T] With K = 0.9 {fQ. ^: For a fixed Af, 
modulation formats optimized for peak power have better 
spectral as well as power efficiency than those optimized for 
average power OOK has the best Po.gain among all studied 
modulation formats. Together with ^, they have a better rj 
than "^^4, 3^p g, ^^^3, and QPSK. The two-dimensional 8- 
ary formats have higher rj than the three-dimensional 16-ary 
formats, and 4-PAM has the highest ry and the lowest Po,gain 
among the studied modulation formats. 

3) Po,gain VS. T] With K = 0.99 (Fig. Q.- As before, 
modulation formats optimized for average power have a higher 
Po,gain than thosc optimized for peak power, ^p 3 and '■^^ 
have the highest Po,gain among the studied formats. QPSK 
and 4-PAM have similar ry which is the highest among the 2-, 
3-, and 4-ary constellations. Unlike the case where K = 0.9, 
the three-dimensional 8-ary constellations are better in spectral 
and power efficiency in comparison to the two-dimensional 8- 
ary constellations. In addition, '^^p ^^ has the highest r] among 
the studied constellations. The two-dimensional constellations 
have more power close to DC, whereas the three-dimensional 
constellations have a wider main lobe. This makes the 99% 
in-band power for the three-dimensional constellations occur 
at frequencies lower than those for the two-dimensional ones. 

4) Po.gain VS. Tj With K — 0.99 (Fig. ^: As in the case 
where K = 0.9, OOK has the highest Po.gain- "^4 has the 
same power efficiency as ,%^, but has a higher 77. QPSK 
and 4-PAM, as in the previous case, have similar ry which 
is the highest among the 2-, 3-, and 4-ary constellations. The 
three-dimensional 8- and 16-ary constellations have higher rj 
than the two-dimensional formats, and the 16-ary constellation 
optimized for peak power has the highest 77. 

V. Conclusions 

We presented a two-dimensional signal space for IM/DD 
that provides a good trade-off between implementation com- 
plexity and spectral efficiency. This signal space suggests sim- 
pler modulator and demodulator structures than for the three- 
dimensional raised-QAM signal space. For a fractional power 
bandwidth of K — 0.9, the two-dimensional formats have 
better spectral characteristics than the three-dimensional ones. 
Therefore, the two-dimensional formats are a good choice for 
single-wavelength optical systems. However, for a fractional 
power bandwidth of K = 0.99, the three-dimensional formats 
have better spectral characteristics. Therefore, they are suitable 
for wavelength-division multiplexing (WDM) systems where 
crosstalk between adjacent channels is important. For both 
signal spaces, modulation formats optimized for peak power 
achieve a higher spectral efficiency than those optimized for 
average power 



Appendix A 
Obtained Constellations 

The coordinates of the optimized constellations in Sec. lIII-Al 
are listed below. The constellations are normalized to have unit 
minimum distance. 

%.3 = {(0, 0), (v^, l/\/3), (^273, -l/Vl)). 



^. 






{(0,0),(V3/2,l/2),(V3/2,-l/2)}. 

X,3U {(2^/273, 0)}. 



^P^^g = ^4 U {(2^273, 2/V3), (2v/273, -2/V3), 
(^6,1/73), (A -l/%/3)}. 



^p„,g = =^4U{(. 



V3/2,l/2 + l/\/3), 



{^/2J?, + %/3/2, -1/2 - 1/V3), 

(2 V273 + V3/2, 1/2), (2v/273 + \/3/2, -1/2)}. 
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